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Abstract.  Er:glass  lasers  have  been  in  operation  with  both  long  pulses 
(hundreds  of  microseconds)  and  Q-switched  pulses  (50  to  100  ns) 
for  more  than  35  yr.  The  ocular  hazards  of  this  laser  were  reported 
early,  and  it  was  determined  that  damage  to  the  eye  from  the 
1 .54-/xm  wavelength  occurred  mainly  in  the  cornea  where  light  from 
this  wavelength  is  highly  absorbed.  Research  on  skin  hazards  has 
been  reported  only  in  the  past  few  years  because  of  limited  pulse 
energies  from  these  lasers.  Currently,  however,  with  pulse  energies  in 
the  hundreds  of  joules,  these  lasers  may  be  hazardous  to  the  skin  in 
addition  to  being  eye  hazards.  We  report  our  minimum  visible  lesion 
(MVL)  threshold  measurements  for  two  different  pulse  durations  and 
three  different  spot  sizes  for  the  1.54-/xm  wavelength  using  porcine 
skin  as  an  in  vivo  model.  We  also  compare  our  measurements  to 
results  from  our  model,  based  on  the  heat  transfer  equation  and  the 
rate  process  equation.  Our  MVL-ED50  thresholds  for  the  long  pulse 
(600  /jis)  at  24  h  postexposure  were  measured  to  be  20,  8.1,  and 
7.4  J  cm"2  for  spot  diameters  of  0.7,  1 .0,  and  5  mm,  respectively. 
Q-switched  laser  pulses  of  31  ns  had  lower  ED50  (estimated  dose  for 
a  50%  probability  of  laser-induced  damage)  thresholds  of  6.1  J  cm"2 
for  a  5-mm-diam,  top-hat  spatial  profile  laser  pulse.  ©  2006  Society  of 

Photo-Optical  Instrumentation  Engineers.  [DOI:  10.1117/1.2187987] 

Keywords:  Er:glass;  laser  skin  model:  skin  injury;  temperature  rise;  laser  safety; 
minimum  visible  lesion,  maximum  permissible  exposure;  estimated  dose  for  a  50% 
probability  of  laser-induced  damage. 
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1  Introduction 

The  Er:glass  laser  has  now  been  in  operation1  with  both  long 
pulses  and  (7-switched  pulses  for  more  than  35  yr.  Ocular 
hazards  of  this  laser  were  reported  early,  and  it  was  deter¬ 
mined  that  damage  to  the  eye  occurred  mainly  in  the  cornea," 
where  light  from  this  wavelength  is  highly  absorbed.  Re¬ 
search  on  skin  hazards  has  been  reported  only  in  the  past  few 
years  because  of  the  limited  pulse  energies  of  Eriglass  lasers 
from  both  short  and  long  pulses  generated  by  these  lasers,  and 
they  were  not  considered  extremely  dangerous.  Currently, 
however,  these  lasers  are  being  produced  with  pulse  energies 
of  hundreds  of  joules  in  the  long-pulse  mode  (millisecond 
pulse  duration),  and  many  joules  for  pulses  in  the  low  nano¬ 
second  range.  Medical  devices3'4  are  being  manufactured  for 
dermatological  uses  such  as  removing  benign  tumors  and  skin 
resurfacing.  Military  applications  include  the  handheld  range 
finders  for  the  U.S.  Army  and  the  low-altitude  navigation  and 
targeting  systems  for  training  purposes  for  the  U.S.  Air  Force 
and  Navy.  The  telecommunication  industry  in  their  fiber  op- 
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tics  system  makes  use  of  this  1 ,54-/xm  wavelength  for  pulse 
modulation  and  transmission. 

Currently,  very  few  skin-damage  threshold  measurements 
have  been  reported  in  the  literature  for  1 . 5 4 - /xm - w a v el c n g t h 
exposures.  Reported  values  vary  to  such  a  wide  extent  that 
additional  measurements  are  required  to  support  the  American 
National  Standards  Institute3  standards  revision.  Experimental 
works  to  extend  the  database  with  consistent  parameters  for 
all  safety  standards,  including  ANSI  Z 136. 1-2000,  Interna¬ 
tional  Electrotechnical  Commission  (IEC)  60825-1,  Food  and 
Drug  Administration  (FDA),  ICNRP,  etc.  are  necessary.  Fur¬ 
thermore,  there  is  a  great  need  for  mathematical  models  that 
can  make  accurate  predictions  of  damage  levels  based  on 
threshold  endpoints  [estimated  doses  for  50%  probability  of 
laser-induced  damage  (ED50s)]  and  reduce  the  required  num¬ 
ber  of  animal  experiments.  It  is  anticipated  that  such  a  model 
would  provide  reliable  predictions  for  any  type  of  laser  for 
different  pulse  durations,  wavelengths,  and  spot  sizes. 

The  ANSI  Z136. 1-2000  standard  divides  the  electromag¬ 
netic  spectrum  into  bands  of  wavelengths.  Wavelengths  from 
1.4  to  100  /rm  are  considered  the  far  IR,  while  the  near  IR  is 
defined  as  the  region  ranging  from  0.7  to  1.4  /jl m.  Maximum 
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permissible  exposures  (MPEs)  are  defined  for  both  the  eye 
and  the  skin  for  the  1 ,54-/xm  wavelength,  and  are  dependent 
on  the  pulse  duration.  This  wavelength  is  referred  to  as  “eye 
safe”  and  has  maintained  that  misnomer  in  the  commercial 
arena  from  its  beginning.  “Eye  safe”  refers  only  to  the  retina 
and  not  the  cornea  or  lens  of  the  human  eye  in  this  context. 
The  MPE  in  the  ANSI  standard  is  set  to  ensure  that  no  injury 
to  the  retina,  cornea,  or  skin  occurs  for  exposures  at  or  below 
the  MPE.  For  example,  the  MPE  for  skin  for  a  1 .54-/xm 
wavelength  exposure  with  a  pulse  duration  ranging  from 
10-9  s  (1  ns)  to  10  s  is  set  at  1  J  cm-2.  Therefore,  for  expo¬ 
sures  of  skin  and  cornea  above  this  MPE  at  1.54  /x m,  severe 
hazards  can  exist,  negating  the  term  “eye  safe.”  For  a  thor¬ 
ough  understanding  of  laser  hazard  evaluation  methods,  see 
Marshall  et  al.6 

The  Yucatan  mini-pig  ( Sus  scrofa  domestica)  was  used  in 
this  study  as  the  model  for  determining  the  threshold  dose  at 
the  minimum  visible  lesion  (MVL)  level  for  laser  injury  to 
skin  at  1.54-/rm  wavelength  for  a  single  0.60-ms  pulse  and  a 
single  31 -ns  pulse.  The  ED50  along  with  the  fiducial  limits  at 
the  95%  confidence  level  were  determined  for  this  exposure 
duration  using  the  Probit7  analysis  technique.  We  used  the 
Yucatan  mini-pig  because  it  has  been  found  to  have  a  higher 
anatomical  similarity  to  human  skin  than  the  commonly  used 
Yorkshire  mini-pig.8  Yucatan  mini-pig  skin  is  melanated  and, 
on  the  flank,  has  a  thickness  similar  to  that  of  the  human  arm, 
which  has  a  high  probability  of  accidental  exposure.  The  data 
obtained  from  this  study  on  porcine  skin  damage  will  contrib¬ 
ute  to  further  understanding  laser  injury  mechanisms.  By  us¬ 
ing  this  model,  the  properties  of  human  skin  can  be  approxi¬ 
mated  more  closely  to  gain  a  better  understanding  of  the 
human/laser-tissue  interaction  for  this  wavelength.  These  re¬ 
sults  will  add  to  the  existing  data  on  laser-skin  effects  on 
which  safety  standards  are  based  and  that  affect  employment 
of  these  lasers. 

2  Materials  and  Methods 

2.1  Micosecond  Regime 

Laser  exposures  were  accomplished  with  a  laser  system 
(Laser  Sight  Technologies  Inc.,  Orlando,  Florida)  using  an 
Er:glass  rod,  delivering  1 ,54-/xm  light  at  a  0.6-ms  exposure 
duration  and  at  various  pulse  energies.  The  laser  produced  a 
Gaussian  beam  profile,  and  it  was  used  for  the  two  smaller 
spot  sizes  (0.7  and  1.0  mm,  measured  at  1/e  irradiance)  in 
this  study.  A  5-mm-diam  spot  size  (top-hat  spatial  profile) 
was  produced  by  a  second  Er:glass  laser  (Megawatt  Lasers, 
75  J/pulse).  Spot  sizes  were  measured  using  an  Electro  Phys¬ 
ics  IR  camera  with  a  Spiricon  LBA  500  laser  beam  analyzer 
and  frame-grabber  card.  The  pulse  duration  was  measured  by 
a  model  ET-3000  InGaAs  Electro-Optics  Technology,  Inc., 
photodiode  connected  to  a  Tektronix  model  TDS  220  oscillo¬ 
scope.  Energy  measurements  were  made  with  a  calibrated 
Molectron  model  JD1000  energy  meter  and  calibrated  J25  and 
150  energy  probes,  which  were  placed  after  a  90/10  beam¬ 
splitter  (BS)  to  collect  10%  of  the  beam  energy,  and  thus 
determine  the  actual  energy  delivered  to  the  skin.  An  articu¬ 
lating  arm  in  the  laser  beam  delivery  system  (Laser  Mecha¬ 
nisms,  Inc.)  was  used  to  deliver  the  beam  without  having  to 
move  the  subject.  A  metal  “aiming  ring”  was  attached  to  the 
end  of  the  articulating  arm,  which  maintained  a  constant  dis- 


Energy 

Probe  90/10  B.S. 


Fig.  1  Block  diagram  of  laser  system  for  delivering  0.6-ms  pulses  at 
1 .54  fj. m  to  skin. 


tance  between  the  arm’s  aperture  and  the  subject.  This  en¬ 
abled  the  precise  positioning  and  distance  control  necessary  to 
deliver  exposures  of  known  spot  size,  more  accurate  beam 
delivery,  and  a  higher  number  of  exposures  per  subject,  result¬ 
ing  in  a  reduction  of  the  total  number  of  subjects  required. 
The  laser  system  setup  is  depicted  in  Fig.  1 . 

2.2  Nanosecond  Regime 

The  Megawatt  laser  used  to  produce  microsecond  pulses  was 
modified  by  installing  an  optomechanical  (/-switch  (Taboada 
Research  Instruments,  Inc.,  San  Antonio,  Texas).  The  modi¬ 
fied  laser  produced  pulses  in  the  30  to  40-ns  range  at  various 
pulse  energies  up  to  3.5  J/pulse.  Pulse  duration  was  measured 
in  the  same  way  as  was  done  in  the  microsecond  regime 
(InGaAs  photodiode  and  oscilloscope).  The  beam  diameter 
for  these  exposures  was  the  same  as  for  this  laser  in  the  mi¬ 
crosecond  regime  (5  mm)  and  was  measured  using  Zap-It 
paper.  The  optical  system  layout  consisted  of  a  90/10  BS 
followed  by  a  1-m-focal-length  lens.  A  final  turning  mirror 
was  used  to  aim  the  beam  to  the  desired  location.  Energy 
measurements  were  made  at  the  location  of  the  BS  using  an 
Ophir  Laserstar  energy  meter  with  Ophir  model  30(150)- 
A-HE  energy  probes.  A  HeNe  laser  was  used  as  marker  beam 
to  locate  the  exposure  point  (see  Fig.  2).  Due  to  the  limited 
spatial  range  of  the  final  turning  mirror,  the  position  of  the 
subject  needed  adjusting  for  some  of  the  exposure  locations. 
A  Panepinto  Sling  (Britz-Heidbrink,  Inc.)  was  utilized  to 
make  the  small  spatial  adjustments.  The  subject’s  elevation 
could  be  changed  as  well  as  the  lateral  position  of  the  subject. 
The  distance  of  the  subject  from  the  laser  was  kept  constant  to 
ensure  uniform  spot  size  for  each  exposure. 

Eleven  female  Yucatan  mini-pigs  (Lonestar  Laboratory 
Swine,  Seguin,  Texas,  and  S  &  S  Farms,  Ramona,  California), 
weighing  between  15  and  20  kg,  were  involved  in  this  study. 
At  least  three  separate  flanks  were  used  at  each  spot  size. 
Animals  were  between  3  and  8  months  of  age.  The  study 
progressed  under  the  animal  use  protocol  titled  “Evaluation  of 
Laser  induced  Corneal  Lesions  in  the  Dutch  Belted  Rabbit 
and  Skin  Lesions  in  the  Yucatan  Mini-Pig,”  which  was  ap¬ 
proved  by  the  Brooks  City-Base,  Texas,  Institutional  Animal 
Care  and  Use  Committee  (IACUC).  None  of  the  animals  were 
euthanized  after  exposure  or  biopsy,  since  they  were  part  of  an 
animal-sharing  program.  Pigs  were  fed  standard  commercially 
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Turning  Mirror  Turning  Mirror 


Fig.  2  Block  diagram  of  laser  system  for  delivering  Q-switched,  31 -ns  pulses  at  1 .54  fj. m  to  skin. 


available  diets  and  had  unlimited  access  to  water.  However, 
solid  food  was  withheld  for  12  h  prior  to  laser  exposure  and 
biopsy  collection.  The  animals  involved  in  this  study  were 
procured,  maintained,  and  used  in  accordance  with  the  Fed¬ 
eral  Animal  Welfare  Act  and  the  “Guide  for  the  Care  and  Use 
of  Laboratory  Animals”  prepared  by  the  Institute  of  Labora¬ 
tory  Animal  Resources,  National  Research  Council.9  Brooks 
City-Base,  Texas,  has  been  fully  accredited  by  the  Association 
for  Assessment  and  Accreditation  of  Laboratory  Animal  Care, 
International  (AAALAC)  since  1967. 

The  Yucatan  mini-pig  subjects  were  received  from  an  at¬ 
tending  veterinary  technician  via  stretcher,  sedated  with  an 
indwelling  intravenous  catheter  placed  in  an  ear  vein,  and 
intubated  prior  to  arrival  at  the  laboratory.  The  pigs  were  se¬ 
dated  by  single  syringe  injection  of  tiletamine/zolazepam  (4  to 
6  mg/kg)  intramuscular  (IM)  and  xylazine  (2.2  mg/kg)  IM, 
and  maintained  on  inhalation  isoflurane  anesthesia  during  all 
procedures.  After  sedation,  hair  on  the  flank  was  clipped  using 
hand  clippers  and  the  underlying  skin  was  cleansed  with  chlo- 
rhexadine  solution  and  allowed  to  air  dry.  The  cleansed  skin 
was  inspected  by  each  of  three  evaluators  to  check  for  red¬ 
ness,  irritation,  or  other  confounding  marks.  All  subjects  were 
kept  warm  during  the  entirety  of  the  procedures.  Physiologi¬ 
cal  parameters  were  monitored  throughout  all  procedures.  Bu- 
prenorphine  (0.05  to  0.1  mg/kg)  was  administered  intramus¬ 
cularly  for  analgesia  after  biopsies  were  complete.  The 
animals  were  returned  to  their  runs  on  recovery  to  sternal 
recumbency  from  anesthesia. 

For  each  subject,  the  flank  to  be  exposed  was  marked  with 
two  6-  X  6-cm  grids  with  a  permanent-ink  marker,  making  a 
total  of  72  grid  squares  per  flank.  As  previously  mentioned, 
the  distance  between  the  articulating  arm  aperture  and  the 
skin  was  kept  constant  by  the  use  of  a  metal  ring  attached  to 
the  end  of  the  arm.  The  animal,  positioned  on  a  table,  did  not 
have  to  be  moved  during  procedures.  Energy  was  delivered 
with  randomly  varying  levels  to  each  square  of  one  grid  and 
this  template  of  energies  was  repeated  on  the  second  grid. 

Postexposure,  three  5 -mm  punch  biopsies  were  obtained 
from  each  subject  immediately  after  the  1-h  postexposure 
reading  and  two  biopsies  again  after  the  24-h  reading.  One  of 
the  three  biopsy  sites  on  each  animal  was  chosen  as  a  control 
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(taken  from  a  location  superior  to  any  exposure  sites).  All 
biopsy  sites  were  closed  with  nonabsorbable  2-0  sutures  and 
topically  medicated  with  Trio-Mycin  ointment  for  infection 
prophylaxis.  Probit7  analysis  was  the  statistical  method  used 
to  determine  the  estimated  dose  for  50%  probability  of  laser- 
induced  damage  (ED50)  for  the  in  vivo  skin  model.  Readings 
of  injured  sites  were  performed  at  1  and  at  24  h  postexposure. 
All  data  points  were  entered  into  the  Probit  statistical  analysis 
package  and  the  ED50s  were  calculated  along  with  their  fidu¬ 
cial  limits  at  the  95%  confidence  level,  slopes,  and  probabili¬ 
ties.  Three  experienced  readers  examined  the  exposure  sites 
and  independently  rated  each  site  for  visible  damage.  At  least 
two  out  of  three  readers  were  required  to  identify  an  exposure 
as  positive  before  a  lesion  was  determined  to  be  present. 

3  Experimental  Results 

Threshold  measurements  for  three  spot  sizes  using  the 
0.6-ms  pulse  duration  are  reported  here.  Enough  data  points 
were  taken  to  provide  the  EDsos  using  Probit  analysis  with 
sufficiently  narrow  fiducial  limits  at  the  95%  confidence  level. 
For  both  the  1-  and  24-h  postexposure  readings,  the  0.6-ms 
ED50s  are  listed  in  Table  1  along  with  their  fiducial  limits  and 
slopes  of  the  Probit  curves  (slope=  Sp/Sd,  where  <5/?= delta 
probability  and  Sd=  delta  dose).  For  the  Q-switched  pulses  of 
3 1  ns,  the  MVL  threshold  was  determined  for  one  spot  size  of 
5  mm  because,  for  the  smaller  beam  diameters,  laser-induced 
breakdown  and  plasma  shielding  occurred  and  prevented  all 
of  the  pulse  energy  from  reaching  the  skin.  Table  2  lists  the 
Q-switched  ED50s  at  1-  and  24-h  postexposure  readings 
along  with  their  fiducial  limits. 

Most  lesions  initially  appeared  as  a  blotched  red  coloring 
near  the  center  of  the  exposure  square,  where  the  laser  beam 
penetrated  the  skin.  These  red  spots  appeared  almost  immedi¬ 
ately  and  many  disappeared  before  the  1-h  reading.  Exposure 
sites  were  observed  after  1  h  and  most  of  the  immediate  le¬ 
sions  were  no  longer  red  splotches  but  very  small  discolora¬ 
tions  in  the  skin.  At  the  24-h  postexposure  reading,  more 
lesions  could  be  clearly  observed  than  were  visible  1-h  pos¬ 
texposure  for  the  two  smaller  spot  sizes.  However,  for  the 
5 -mm  spot  size,  fewer  lesions  were  observed  at  the  24-h 
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Table  1  MVL-ED50  for  a  0.6-ms  pulse  duration  of  1  .54-/rm  laser  pulses. 


Experimental  Setup, 

Number  of  Subjects,  and  Shots 

MV1-ED50  (J  cm"2) 

1  -h  Reading 

MVL-ED5o  (J  cm"2) 
24-h  Reading 

Probit  Curve 
Slope=  Sp/ Sd 

0.7-mm-diam  spot 

4  pigs,  4  flanks, 

263  exposures 

26  (30-23) 

20  (21-18) 

5.7 

1 .0-mm-diam  spot 

3  pigs,  3  flanks, 

2 1  6  exposures 

1  1  (1  1-9.5) 

8.1  (8. 7-7.5) 

15 

5. 0-mm-diam  spot 

4  pigs,  4  flanks, 

287  exposures 

6.4  (6. 7-6. 2) 

7.4  (7. 8-7.0) 

14 

reading  and  thereby  producing  a  larger  ED50  at  24  h  than  at 
the  1-h  reading.  For  this  spot  size  only,  the  threshold  was 
lower  and  the  fiducial  limits  had  a  reduced  spread. 

Additional  Probit  runs  were  made  for  the  0.7-mm-diam 
spot  size  when  it  was  noticed  that  the  probability  of  chi- 
square  for  the  24-h  reading  was  a  very  low  value  of  0.0057 
while  the  Pearson’s  chi-square  was  a  value  of  208.  For  the 
279  exposures  read  at  24  h.  the  radiant  exposures  ranged  from 
4.9  to  5 1  J  cm-2  and  it  was  determined  that  several  zeros  (no 
visible  lesions)  were  recorded  for  the  high-energy  exposures. 
A  subset  of  data  was  run  with  all  exposures  of  42  J  cm-2  and 
above  removed  from  the  table,  which  gave  the  same  threshold 
but  with  the  chi-square  increased  to  0.48.  This  same  subset  of 
data  was  again  revised  to  limit  the  maximum  to  33  J  cm-2  for 
the  top  exposure,  and  the  chi-square  was  again  increased  to 
0.89  without  changing  the  threshold  of  20  J  cm-2.  Pearson’s 
chi-square  was  decreased  from  208  to  108  for  this  run.  The 
ED50  values  determined  for  each  of  these  sets  did  not  vary  by 
more  than  10%,  and  the  fiducial  limits  overlapped  for  all  three 
conditions.  These  changes,  when  applied  to  the  1-h  reading, 
did  not  significantly  affect  any  parameters.  Since  the  chi- 
square  for  the  1-  and  5-mm  beam  diameters  were  very  close 
to  1.00  at  the  1-  and  24-h  reading,  they  are  reported  as  origi¬ 
nally  run  with  all  the  exposures.  These  results  are  plotted  in 
Fig.  3  together  with  other  referenced  data. 

4  Thermal  Modeling 

Our  mathematical  model  (Takata  Thermal  model)  was 
developed10'11  at  Brooks  Air  Force  Base  (AFB),  Texas,  in  the 
early  1970s.  The  model  is  currently  being  run  on  commodity 
personal  computers  running  a  Linux  operating  system,  using 
legacy  FORTRAN  source  codes.  The  Takata  Thermal  Model 
consists  of  three  separate  validated  capabilities:  one  model  for 
the  retina,  a  second  for  the  cornea,  and  a  third  model  for  skin. 


All  three  implement  a  finite-difference  numerical  solution  of  a 
time-dependent  2-D  bioheat  equation  in  cylindrical  coordi¬ 
nates,  following  an  early  work  by  Mainster  et  al.12  on  heat 
transfer  in  biological  media.  The  result  is  a  simulation  of  the 
temperature  rise  as  a  function  of  time  at  each  point  in  the 
tissue.  The  models  are  coupled  to  a  rate-processes  damage 
integral  [Eq.  (1)]  to  obtain  an  evaluation  of  potential  damage 
to  the  tissue.  The  models  were  originally  developed  for  the 
visible  wavelength  bands  (for  nonhuman  primate  eyes).  They 
were  later  modified  and  validated  for  nonporcine  skin11'13  and 
nonhuman  corneas  through  the  implementation  of  boundary 
conditions  appropriate  for  an  air-tissue  interface.  They  have 
since  been  revised  and  extended  to  include  various  IR  wave¬ 
lengths  and  temperature  effects  such  as  charring  of  tissues. 
Calculations  using  this  model  can  be  performed  on  all  three 
tissue  types  for  virtually  any  wavelength  for  which  optical 
properties  of  the  tissues  are  known.  One  of  the  most  challeng¬ 
ing  issues  in  the  use  of  these  models  for  the  near-IR  region  is 
the  paucity  of  key  parameters  beyond  1200  nm  (i.e.,  absorp¬ 
tion  and  scattering  coefficients). 

For  this  study,  skin  modeling  capability  was  employed  to 
estimate  temperature  effects  along  with  the  evaluation  of  tem¬ 
perature  rise  distributions  and  thresholds  for  tissue  damage. 
Features  of  the  model  include  a  user-configurable  multilayer 
tissue  model.  Thermomechanical  as  well  as  optical  properties 
of  the  skin  tissue  layers  are  user  inputs.  Laser  parameters  are 
also  user  configurable  and  a  geometric  model  of  beam  irradi- 
ance  is  employed  along  with  linear  absorption  of  the  tissue  to 
estimate  energy  deposition  rates  at  various  points  within  the 
computational  grid.  Boundary  conditions  include  constant 
flux  surface  convection  at  the  tissue-air  interface.  Thermal 
effects  of  variable  blood  flow  with  tissue  depth  are  evaluated. 
The  phase  change  of  the  water  content  of  the  tissue  as  well  as 
increased  absorption  for  charred  tissues  are  evaluated  through 


Table  2  MVL-ED50  data  for  1 .54-/nn  laser  pulses  at  a  Q-switched  pulse  duration  of  31  ns. 


Experimental  Setup, 

Number  of  Subjects,  and  Shots 

MV1-ED50  (J  cm"2) 

1  -h  Reading 

MVI-ED50  (J  cm"2) 
24-h  Reading 

Probit  Curve 
Slope  =  Sp/  Sd 

5. 0-mm-diam  spot 

6.3  (6. 8-5. 8) 

6.1  (6. 5-5. 5) 

6.6 

2  pigs,  3  flanks, 

2 1  6  exposures 
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Model  Threshold  (for  2nd  &  3rd  Degree  Burn)  Versus  Spot  Size 


Fig.  3  Thresholds  for  second-  and  third-degree  burns  versus  spot  size. 

empirical  relationships.  The  model  does  not  incorporate  tissue 
optical  scattering  effects.  However,  we  estimate  that  for  beam 
diameters  of  less  than  100  /rm,  this  would  be  a  significant 
effect  for  this  short-pulse  study  at  near-IR  wavelengths. 

The  Takata  Thermal  Model  execution  results  in  a  time- 
temperature  history  at  each  point  within  the  computational 
grid.  Each  point  within  the  grid  is  evaluated  for  potential 
damage  over  the  duration  of  the  simulation  through  an 
Arrhenius  damage  integral  with  temperature  dependent  dam¬ 
age  rate  coefficients. 

The  damage  model12  we  used  is  based  on  the  work  of 
Henriques14  in  pig  skin,  and  takes  the  form  given  by 

D,(r,z)=A  f  exp(- E/RT)  dr,  (1) 

Jo 

where  A  is  a  preexponential  normalization  factor  (in  inverse 
seconds),  E  is  the  activation  energy  for  the  reaction  (cal/M),  R 


is  the  universal  gas  constant  (2.0  cal/M  K),  T=T(r,z,t)  is 
the  locale  absolute  temperature  (in  kelvins),  and  t  is  the  time 
at  final  recovery  of  temperature  after  exposure.  Henriques 
used  rate  coefficients  for  Eq.  (1)  such  that  complete  necrosis 
of  the  skin  was  associated  with  17  =  1 ,  and  determined  that 
A  =  3.1  X  1098  (1/s),  and  E=  150,000  (cal/M).  The  Takata 
Skin  Model  uses  values  that  vary  as  a  function  of 
temperature,11  establishing  multiple  rates  of  damage.  Values 
of  A  =  4.32  X  1064  (1/s)  and  E=  100,000  (cal/M)  for  tem¬ 
perature  ranges  of  317  to  323  K,  and  A  =  9.389 
Xl0104(l/s)  and  E=  160,000  (cal/M)  for  temperature 
ranges  of  323  to  333  K.  Below  317  K,  the  damage  contribu¬ 
tion  is  assumed  to  be  zero. 

In  the  Takata  (Skin  Damage)  Thermal  Model,  the  damage 
integral  is  normalized  against  experimental  data  for  first- 
through  third-degree  burns.  For  the  purposes  of  this  paper,  the 
second-degree  burn  is  represented  by  a  damage  integral  value 
of  0=1.0  and  the  third-degree  burn  by  a  value11  of  17 
=  10,000. 


5  Modeling  Results 

Calculations  from  the  Takata  Thermal  Model11  for  skin  are 
listed  in  Table  3  for  minimal  second-degree  burns  calculated 
for  various  pulse  energies  and  spot  sizes  using  only  the  long 
pulse  (0.6  ms).  Laser  exposure  diameters  were  varied  from 
0.1  to  5  mm  and  the  laser  pulse  powers  were  varied  to  deter¬ 
mine  the  power  required  to  provide  a  minimal  degree  bum. 
Minimum  temperature  rises  are  plotted  as  a  function  of  spot 
diameter  in  Fig.  4  for  those  values  listed  in  Table  3  calculated 
for  this  pulse  duration.  The  radiant  exposure  required  to  pro¬ 
duce  these  temperature  rises  are  listed  in  Table  3  and  are 
plotted  in  Fig.  3  together  with  the  MVL-ED50  data  at  24  h 
listed  in  Table  1  and  two  other  data  points  from  the  literature. 
In  this  figure,  it  can  be  seen  how  closely  the  Takata  Skin 
Model  computes  the  thresholds  as  measured  in  the  laboratory 
for  radiant  exposure  values  with  a  spot  diameter  equal  to  or 
greater  than  1.0  mm. 


Table  3  Takata  Second-Degree  Burn  Model  outputs  for  eight  spot  sizes  using  latest  parameters  for  skin 
at  a  pulse  duration  of  0.6  ms. 


Spot  Diameter 
(mm) 

Pulse  Energy  (J) 

Power  (W) 

Temperature  (°C) 

Radiant  Exposure  (J  cm  2) 

0.10 

0.00126 

2.1 

43.2 

16.04 

0.25 

0.00678 

11.3 

40.1 

13.81 

0.50 

0.0246 

41.0 

36.9 

12.53 

0.70 

0.0462 

77.0 

35.5 

12.00 

1.00 

0.0906 

151 

33.9 

1  1.54 

2.00 

0.3324 

554 

31.1 

10.58 

3.00 

0.72 

1200 

30.0 

10.19 

4.00 

1.26 

2100 

29.5 

10.03 

5.00 

1.95 

3250 

29.2 

9.93 
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Model  Temperature  Rise  Vs.  Spot  Size  For  2nd  Degree  Burn 


Fig.  4  Takata  Model  temperature  rises  for  second-degree  burns  for 
various  spot  diameters. 


We  performed  the  same  calculations  using  the  31 -ns  pulse 
duration  for  a  spot  diameter  of  5  mm  using  the  ED50  value  of 
power  and  the  model  predicted  there  would  be  no  damage. 
The  temperature  rise  calculated  was  only  4.7  °C  with  the 
ED50  pulse  energy  and  the  damage  level  calculated  to  be  zero. 
We  increased  the  pulse  power  enough  to  give  a  second-degree 
burn  for  a  long  pulse  and  the  required  pulse  energy  was  5 
times  the  threshold.  At  this  pulse  energy,  the  temperature  rise 
was  calculated  to  28.6° C.  This  temperature  rise  compares 
well  with  the  29.2° C  rise  calculated  for  a  0.6-ms  pulse  at  the 
ED50  energy  (see  Table  3). 

6  Discussion 

In  our  data  presented  here,  three  different  spot  sizes  for  a 
pulse  duration  of  0.6  ms  at  1.54  /rm  with  ED50  measure¬ 
ments  are  reported  with  a  comparison  to  other  reported  mea¬ 
surements  and  to  a  mathematical  thermal  model.  We  also  re¬ 
port  on  (I-switchcd,  nanosecond  pulses  for  a  spot  size  of 
5  mm  and  make  the  same  comparison  with  reported  values. 
Laser-induced  damage  to  porcine  skin  has  been  reported  by 
Rico  et  al. 15,16  for  the  1 ,54-/xm  wavelength  and  a  small  spot 
size  (0.57-mm  diameter)  and  they  used  two  pig  models 
(Yorkshire  and  the  Yucatan  mini-pigs)  for  their  threshold  mea¬ 
surements.  They  used  a  smaller  spot  size  (0.57  mm)  than  the 
one  used  in  this  study  because  they  were  limited  in  the  avail¬ 
able  pulse  energies  from  their  laser.  They  reported  42  J  cm-2 
for  the  Yorkshire  and  53  J  cm-2  for  the  Yucatan. 

We  compare  other  data  points  with  ours  in  Fig.  3  that 
indicates  that  there  is  a  strong  spot  size  dependency  for  expo¬ 
sure  sites  below  1  mm  in  diameter.  Rico’s  threshold  value 
was  considerably  larger  than  our  values  for  the  Yucatan  mini¬ 
pig  for  spot  sizes  of  0.7  and  1 .0  mm  (20  J  cm-2  for  0.7  mm 
and  8.1  Jem-2  for  1  mm  diameter).  Similar  measurements 
were  reported  by  Lukashev  et  al. 17,18  but  with  larger  spot  sizes 
using  white  “Krupnaya  belaya”  pig  skin.  They  report  the  dam¬ 
age  threshold  for  a  spot  size  of  5.5  mm  in  diameter  for  pulse 
durations  of  2.5  ms  at  24  h  to  be  6.5  J  cm-2  and  a 
2.5-mm-diam  beam  for  100  ns,  giving  3.5  J  cm-2  at  24  h. 
Their  2.5-ms  threshold  6.5  J  cm-2  is  also  shown  in  Fig.  3  for 
a  comparison  with  our  5 -mm  spot  size. 
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For  the  three  thresholds  listed  in  Table  1  at  the  24-h 
postexposure,  two  showed  a  lower  threshold  than  at  the  1-h 
reading.  All  three  thresholds  shown  in  Fig.  3  are  at  the  24-h 
reading  and  two  of  the  three  had  a  lower  threshold  than  pre¬ 
dicted  by  the  Takata  Model  for  a  second-degree  burn.  Natu¬ 
rally,  the  third-degree  bum  thresholds  shown  in  Fig.  3  would 
require  a  higher  fluence  than  the  second-degree  bum  and 
these  data  are  included  only  to  show  relative  levels.  The  other 
two  data  points  (Rico  and  Lukashev)  are  included  to  indicate 
the  range  of  values  to  be  expected  for  large  and  small  spot 
sizes.  Our  data  agree  with  Lukashev’s  data  at  the  larger  spot 
size,  but  not  so  well  with  Rico’s  data  at  a  smaller  spot  size. 
Since  our  threshold  value  for  the  0.7-mm  spot  size  was  so 
much  larger  than  the  model  predictions  and  the  fact  that  we 
were  able  to  see  flashes  of  light  at  the  delivery  site  and  some¬ 
times  hear  a  “pop,”  we  believe  that  breakdown  was  occurring 
at  or  near  the  surface  of  the  skin  and  not  all  of  the  pulse 
energy  was  penetrating  the  skin.  In  analyzing  our  data,  we 
discovered  that  there  were  several  zeros  or  nonlesions  at  these 
very  large  laser  pulse  energies.  In  fact,  the  chi-square  term  in 
the  Probit  analysis  was  found  to  be  0.007,  showing  this  dis¬ 
tribution  deviated  highly  from  a  log-normal  distribution.  We 
determined  it  appropriate  to  remove  all  data  too  far  above 
threshold  or  data  above  a  certain  exposure  level  from  the  Pro¬ 
bit  analysis  because  of  the  likelihood  of  generation  of  plasma 
that  shielded  the  skin  from  the  laser  energy,  thereby  giving  a 
false  negative  reading.  When  we  removed  these  data  points 
above  certain  pulse  energies,  we  raised  the  chi-square  term  to 
above  0.6,  but  the  ED50  value  did  not  change. 

We  hypothesize  that  some  of  these  pulse  energies  were 
much  too  large  for  the  spot  diameter  being  used  and  that  this 
energy  was  probably  creating  plasma  at  the  skin  surface  in¬ 
stead  of  propagating  into  the  skin  where  it  could  contribute  to 
the  creation  of  a  lesion.  Thus,  some  of  the  pulse  energy  could 
have  been  absorbed  before  the  pulse  reached  the  surface  of 
the  skin  and  this  would  have  required  a  higher  energy  pulse  to 
cause  a  visible  lesion.  To  date,  we  have  found  no  other  reason 
for  such  a  large  deviation  between  model  thresholds  and  the 
ED50  values  for  spot  diameters  below  1 .0  mm.  We  note  that 
the  chi-square  distribution  calculated  for  the  1 .0-mm  spot  di¬ 
ameter  was  1.00,  and  flashes  of  light  normally  associated  with 
laser-induced  breakdown  (LIB)  production  were  not  observed 
at  this  spot  diameter.  We  again  hypothesize  that  the  threshold 
reported  by  Rico  et  al..16  also  had  LIBs  and  plasmas  generated 
due  to  its  very  high  threshold  and  conclude  that  more  data 
points  below  1  mm  in  diameter  should  help  to  clarify  this 
ambiguity. 

Our  thermal  model  predicts  a  temperature  rise  due  to  en¬ 
ergy  deposition  within  the  tissue  throughout  a  volume  as  de¬ 
fined  by  the  input  parameters.  In  the  Takata  Skin  Model,  this 
prediction  of  temperature  rise  is  based  on  power  input  where 
the  model  computes  a  time  step  increase  throughout  the  en¬ 
ergy  deposition  (laser  pulse  duration)  and  for  as  long  after¬ 
ward  as  necessary  to  return  the  temperature  back  to  preexpo¬ 
sure  levels.  We  note  that  this  spot-size  dependency  is  greater 
for  spot  sizes  less  than  1 .0  mm  in  diameter  as  seen  in  Figs.  3 
and  4.  Above  a  1.0- mm  spot  diameter,  the  Takata  Skin  Model 
temperature  rises  have  been  correlated  with  damage  levels  as 
a  validation  of  the  thermal  model,  while  below  1 .0  mm,  no 
correlation  is  yet  possible. 
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Fig.  5  Pigskin  showing  Q-switched,  1 .54-/r,m  superthreshold  lesions 
2  weeks  postexposure. 


For  the  Takata  Skin  Model,  Fig.  3  shows  the  data  of  power 
and  spot  size  in  terms  of  the  radiant  exposure  in  joules  per 
square  centimeter  necessary  to  produce  the  minimal  observ¬ 
able  second-degree  burn  injury  on  the  area  of  the  exposure. 
For  spot  sizes  greater  than  1  mm  in  diameter,  there  is  very 
little  dependence  on  the  area  and  the  threshold  irradiance  re¬ 
mains  essentially  constant.  In  contrast  to  actual  ED50s  below 
1 .0  mm,  the  Takata  values  show  a  slight  dependence  on  laser 
spot  size.  The  interdependence  of  time  and  temperature  is 
critical  with  any  rate-process  model.  Therefore  this  depen¬ 
dence  on  the  time-temperature  history  and  not  on  temperature 
rise  alone  must  be  worked  out  for  laser  spot  input  diameters 
less  than  1 .0  mm. 

For  the  nanosecond  pulses,  the  ED50  thresholds  given  in 
Table  2  are  very  nearly  the  same  as  for  the  0.6-ms  pulses  at 
the  same  spot  diameters.  It  was  not  anticipated  that  they 
would  be  the  same,  but  considering  the  results  of  Lukashev  et 
al.  for  the  white  pig,  the  results  are  not  that  different.  They 
reported  a  threshold  of  3.5  J  cm-2  for  the  100-ns  pulse  as 
compared  to  the  6.5  J  cm-2  for  the  millisecond  pulses.  Their 
spot  size  was  only  half  the  size  of  our  5 -mm  diameter  and 
this,  as  well  as  the  differences  in  pigskins,  may  have  contrib¬ 
uted  to  the  lower  threshold.  When  we  tried  to  use  smaller  spot 
sizes,  there  was  breakdown  at  the  skin  surface  and  plasma 
shielding,  giving  a  flash  of  light  and  a  loud  “pop”  sound.  Even 
at  the  5 -mm  beam  diameter,  with  pulse  energies  more  than 
twice  the  ED50,  some  pulses  produced  breakdown. 

After  two  weeks,  most  of  the  superthreshold  lesions  were 
still  visible  on  the  skin.  Figure  5  shows  some  of  these  lesions 
for  different  energy  levels  from  barely  visible  to  large  brown 
spots.  The  large  black  spot  near  the  center  of  the  photo  is  the 
remains  of  the  biopsy  sample  that  shows  the  healing  process. 
At  this  time  histology  studies  are  ongoing  in  a  larger  study 
and  will  not  be  reported  in  this  work. 

When  we  computed  the  temperature  change  using  the  ther¬ 
mal  model,  the  increase  in  temperature  was  only  4.7 °C;  not 
nearly  enough  to  give  the  thermal  damage  predicted  by  the 
rate  equation  used  in  the  model.  This  appears  to  be  reasonable 
when  one  considers  that  the  energy  delivered  to  the  tissue  is 
defined  as  power  multiplied  by  time,  and  as  the  time  becomes 
shorter,  the  power  increases  for  like  energies.  The  small  tem¬ 
perature  increase  observed  in  the  model  was  related  to  very 
high  peak  powers  found  in  nanosecond  pulses,  as  compared  to 
longer  pulse  durations  such  as  microsecond  pulses.  In  the 
nanosecond  time  regime,  it  is  possible  that  steam  bubbles  can 
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be  generated  at  the  surface  of  the  skin,  and  thereby  reduce  the 
energy  penetrating  the  skin  to  cause  damage.  We  recognize 
that  our  model  does  not  predict  a  sufficient  temperature  rise  to 
cause  injury  at  the  experimentally  measure  threshold.  It  is 
obvious  that  a  pure  heat  transfer  model  may  not  be  sufficient 
to  accurately  predict  damage  thresholds  when  an  ionization 
(plasma  formation),  ablation,  or  microcavitation  event  occurs 
within  the  tissue  as  the  threshold  damage  mechanism.  Further 
work  will  be  required  to  incorporate  all  physical  mechanisms 
into  a  combined  physical  model. 

Thus,  we  believe  that  the  damage  observed  was  caused  by 
acoustic  or  shock  waves  rather  than  thermal  denaturation.  It  is 
interesting  to  note  that  we  calculated  the  temperature  rise  for 
the  same  pigskin  using  50-ns  pulses  at  1318  nm  and  found 
similar  results,  only  a  few  degrees  temperature  rise  and  zero 
damage  for  the  ED50  pulse  energies  at  that  wavelength.  Luka¬ 
shev  et  al.  calculated  a  temperature  rise  of  8  °  C  for  nanosec¬ 
ond  pulses  and  15  °C  for  millisecond  pulses  at  their  ED50 
values,  comparing  favorably  with  our  temperature  calcula¬ 
tions.  They  used  an  absorption  coefficient  of  a=10cm-1, 
while  we  used  8  cm-1  in  our  model.16 

We  can  also  compare  the  preceding  thresholds  with  those 
measured  for  ablation  using  nanosecond  laser  pulses  at  other 
wavelengths.  Hu  et  al.19  studied  the  ablation  thresholds  for 
different  wavelengths  from  213  to  1064  nm  for  14-ns  pulses 
and  reported  that  the  thresholds  varied  from  1 .6  J  cm-2  at 
213  nm  to  129  J  cm-2  at  1064  nm.  These  values  were  for  the 
EDqo  thresholds  and  the  closest  wavelength  to  our  1540  nm 
was  1064  nm.  They  reported  using  an  absorption  coefficient 
a  for  white  domestic  pigskin  of  5  cm-1  compared  to  our 
value  of  a  of  8  cm-1  for  the  Yucatan  mini-pig.  Using  prob¬ 
ability  values  of  EDg0  to  compare  with  their  thresholds,  we 
calculated  a  9.5  J  cm-2  value,  which  is  15  times  smaller  for 
the  visible  lesion  ED50,  compared  with  the  ablation  threshold 
for  a  I  -fim  ablation  depth.  They  focused  the  beam  to  a  spot 
diameter  of  14  /rm  as  compared  to  our  5-mm-diam  spot  size. 
Thus,  our  threshold  would  have  been  much  greater  had  we 
used  such  a  small  laser  beam  diameter. 

The  MPE  as  set  by  ANSI  for  skin  exposures  is  a  constant 
value  of  1  J  cm-2  for  pulse  durations  from  1  ns  to  10  s,  and 
this  value  covers  both  of  the  pulse  duration  reported  herein. 
For  our  larger  spot  sizes  of  5-mm  diameter,  we  measured  7.4 
and  6.1  J  cm-2,  which  are  both  greater  than  the  MPE.  Theses 
values  are  comparable  to  values  of  6.5  and  3.5  J  cm-2  re¬ 
ported  by  Lukashev  for  the  millisecond  and  nanosecond  pulse 
durations,  respectively. 

7  Conclusions 

We  reported  data  that  may  be  included  in  the  databank  for 
skin  threshold  measurements  and  these  data  are  compared 
with  the  ANSI  Z136. 1-2000  standard  for  skin  MPEs  at 
1.54  /rm  and  other  reported  skin  measurements.  We  also 
compare  these  to  thermal  damage  model  calculations  for 
second-degree  burns  to  the  skin  for  laser  spot  sizes  from 
0.1  to  5  mm  in  diameter.  We  used  the  standard  finite- 
difference  thermal  model  of  heat  conduction  (Takata  et  al.) 
coupled  with  Henriques’s  rate  process  model  of  thermal  dam¬ 
age.  With  this  model,  we  can  determine  the  maximum  tem¬ 
perature  reached  during  and  after  the  laser  pulse  for  pulses 
longer  than  nanoseconds,  temperature  distribution  throughout 
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the  tissue  and  the  damage  levels  reached  within  tissue.  The 
“American  National  Standard  for  Safe  Use  of  Lasers”  (ANSI- 
Z136.1,  2000)  gives  the  skin  MPE  for  the  far-IR  wavelengths 
(1.5  to  1.8  /rm)  as  1.0  Jem-2  for  pulse  durations  of 
1  ns  to  10  s.  For  spot  sizes  less  than  1  mm  in  diameter,  our 
ED50  for  the  Yucatan  mini-pig  skin  was  20  times  the  MPE  and 
for  the  1  -mm  diameter,  it  was  still  8  times  the  MPE.  For  spot 
sizes  greater  than  1  mm,  our  data  show  that  the  threshold  was 
still  7.4  times  the  MPE,  and  this  compares  favorably  with  the 
other  reported  data  point  for  a  5.5-mm  spot  of  6.5  J  cm-2  by 
Lukashev  et  al.  For  the  nanosecond  pulses,  our  MVL  thresh¬ 
old  was  6  times  the  MPE  and  it  was  15  times  below  the 
reported  ablation  thresholds  (Hu  et  al.)  for  near-IR  wave¬ 
lengths.  Thus,  we  conclude  that  the  MPE  as  set  by  ANSI  is 
below  all  of  our  measurements  and  probably  does  not  require 
revision  unless  new  and  different  data  are  reported. 
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